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Using the density matrix renormalization group algo- 
rithm (DMRG), we study the model of spinless fermions with 
nearest-neighbor interaction on a ring in the presence of dis- 
order. We determine the spatial decay of the density induced 
by a defect (Friedel oscillations), and the phase sensitivity of 
the ground state energy aE = { — )^{E{(j) = 0) — E{(j) = tt)), 
where — 27r"l>/<E>o {N is the number of fermions, $ the mag- 
netic flux, and "l>o = h/e the flux quantum), for a disordered 
system versus the system size M. The quantity In (A'IaE) is 
found to have a normal distribution to a good approximation. 
The "localization length" decreases (increases) for a repulsive 
(attractive) interaction. 



INTRODUCTION 

About five decades after the first theoretical discussion 
of orbital magnetism of free electrons on a ring , this 
phenomenon — or "persistent" current - has been ob- 
served in mesoscopic metallic and semiconducting 
ring structures. At present, the magnitude of the 
effect — experimental results are much larger than the- 
ory predicts — is not well understood. It was suggested 
immediately after the experiment that the constraint of 
local charge neutrality, imposed by the electron-electron 
interaction, could lead to a considerable enhancement of 
the current [^,0, but these results, at least for single 
rings, were not conclusive B, and were not confirmed 
numerically ^ (see also e.g. [ |lO|jTl| ]). 

Considering electrons on a ring, the ground state en- 
ergy, E{(j)), depends on the boundary condition, charac- 
terized by the phase (0 = corresponds to periodic, 
and (j) — ±7r to antiperiodic boundary conditions). Al- 
ternatively, the boundary condition can be interpreted as 
arising from a magnetic flux, <i>, provided we identify (j) 
with 27r<l>/$o, where = h/e denotes the flux quantum. 
Clearly, E{(j)) is periodic with period 27r. The energy dif- 
ference between periodic and antiperiodic boundary con- 
ditions, aE, the persistent current, / ^ —E'{(j)), and the 
charge stiffness (the "Drude weight"), D - E"{(j) = 0), 
are a measure of the phase sensitivity of the system (a 
comprehensive discussion is given in |]l2| ) . In view of the 
discrepancy between experiment and theory, in particular 
for the metallic samples [|]J^, it is important to under- 
stand further the interplay between interaction and disor- 
der. Progress has been made for simple one-dimensional 
models on these we concentrate in the following. 

We investigate the standard model of spinless fermions 



on a ring, with nearest-neighbor interaction V (in units 
of the hopping amplitude), the lattice constant is unity, 
M denotes the number of sites (i.e. the system size), and 
N the fermion number. We restrict ourselves to the case 
of half-filling, N = M/2, and consider on-site disorder. 
The model is described by the Hamiltonian H = Hk + 
Hv + 'Hi, with 



where = (p/M , 



and 
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In contrast to the "generic" defects described by the 
random on-site energies {ej}, it is also possible to con- 
struct integrable models with "transparent" impu- 
rities, which decrease the phase sensitivity even though 
there is no wave-function localization. The distinction 
between "integrable" and "non-integrable" is also most 
important for the temperature dependence of the Drude 
weight, D We remark that the clean case is well 

studied by exact methods (see e.g. |la,M). 



THE DMRG-ALGORITHM 

The Density Matrix Renormalization Group algo- 
rithm, introduced by White |jl^, is a numerical tech- 
nique which allows reliable results for one-dimensional 
quantum lattice models, of a size (up to a few hundreds 
sites) much larger than accessible by exact diagonaliza- 
tion methods, to be obtained. The algorithm can, in 
short, be characterized as a "projected diagonalization" , 
where the subspace to be projected onto is determined 
by the most probable eigenstates of a density matrix. 
Start, for example, with a reasonable representation of a 
s-site system, using m relevant states. Then add one site, 
s — > s -|- 1, and supplement the system by an "environ- 
ment" , namely the sites s + 2, . . . , 2(s + 1). The basis of 
(2rn)^ states formed in this way is used to determine the 
ground state of H2s+2- (The factor "two" appears here 
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since we have two states for each site.) The density ma- 
trix of the "system" determines the m most important 
states, onto which all relevant operators are projected. 
Then add another site, s + 1 s + 2, and proceed. 

The result of this "infinite lattice" algorithm is the ba- 
sis for the "finite lattice" algorithm, which is needed in 
order to treat non-reflection symmetric systems, (i) The 
infinite lattice algorithm is used to find a representation 
of all s-sites systems with s up to M — 3. (ii) Then con- 
sider the sites 1, . . . , s as system, add two sites, and take 
as environment the s + 3, . . . , M system, (iii) Proceed 
as described above. This process can be characterized as 
"sweeping through the lattice". With m up to several 
hundred, we have achieved sufficient accuracy, even for 
questions as subtle as the phase sensitivity. The method 
works best for open boundary conditions, though there is 
no general problem (except enhanced computing time) to 
include twisted boundary conditions. Since the DMRG 
is a local method, disorder is easily included, though dis- 
order averages of course require considerable computing 
time. For example, the data given in Fig. ^ are based 
on roughly 300 CPU days on a high-end workstation. 
Further details, e.g. the adaption to non-reflection sym- 
metric models, are discussed in [H. 



FRIEDEL OSCILLATIONS 

The decay of the density oscillations induced by a de- 
fect is a long-standing problem in solid state physics. 
This phenomenon, called Friedel or Ruderman-Kittel os- 
cillations (depending on the context), is closely related to 
the singularity in the response function for wave-vectors 
close to 2kF- It is expected that asymptotically, the in- 
duced density decays as 



5n(x) 



cos {2kpx + rjp) 



(4) 



Using the DMRG, we have computed 5n{x) for a sys- 
tem of 200 sites and various interaction strengths and, 
as a test for the accuracy of our calculation, for systems 
with M = 500 and V = ±f,2. The impurity is chosen 
antisymmetrically for technical reasons, ei = —cm- We 
consider a half- filled band, i.e. 2kpXj = nj, where j is 
the distance (in units of the lattice spacing) from the de- 
fect. A sample of our results is shown in Fig. 0, where we 
plot, on a logarithmic scale, the magnitude \Snj\ versus 
distance. Clearly, it is possible to extract the exponent 
S without difficulty. We emphasize that the algebraic 
decay starts already at a few lattice site sites. Typi- 
cally, we have used a basis of ni — 120 (200) states for 
the M = 200 (500) system, and performed four sweeps 
through the lattice. 

The exponent 6 as a function of the nearest-neighbor 
interaction V is given in Fig. |^. The exponent S decreases 
with increasing repulsive interaction, and increases with 



attractive interaction, compared to the value for non- 
interacting fermions, 6 = 1 (one dimension!). 

Qualitatively, this trend agrees with the predictio n |19| ] 
based on the Luttinger liquid. In a recent preprint [|20{, S 
was related to the "dressed charge" of the (clean) model, 
with the result S = = 7i'/4?7, where rj, related to V 
through V = —2cos{2ri), parameterizes the interaction. 
The expression 6 = Tr/Arj is also shown in Fig. ^ and 
is in almost perfect agreement with our numerical data, 
except for y > 1 where we find that the oscillations decay 
more weakly than predicted. This seems to be related to 
the crossover (for a weak impurity, and V > 0) found in 
Ip^ , i.e. for the system sizes studied we may not yet be 
in the asymptotic regime. We have preliminary results 
showing that for a strong impurity, S tends to increase 
towards the asymptotic result given in [E0|. 



PHASE SENSITIVITY 

Concerning the phase sensitivity of the ground state 
energy, let us recall the free electron result, for odd N: 



[i?(0)-i?(O)]o = g|0^ 



(5) 



to be continued periodically outside the interval — tt . . . tt. 
For an even number of particles, —><?!) — tt in this equa- 
tion. In our units, and for a half-filled band, we have 
hvp/L 2/M- thus aEq = (-1)^[£;(0) - E{it)]v=() = 
■K /M for the clean, non-interacting system. 

Disorder is introduced by taking the on-site energies 
{ej } as random quantities, uniformly distributed over the 
range —14^/2 . . . 14^/2, which corresponds for free fermions 
to a conductance localization length ~ 105/ Vt^^, i.e. 
the average conductance decreases {g) ^ exp(— 2Af/^o) 
for a wire of length M |^ . We considered W = 2 only, 
hence ~ 26. Naturally, we first studied non- interacting 
fermions, using 500 samples, i.e. 500 different realizations 
of the disorder for each system size. Interestingly, we find 
(aEq) ^ (g)^^^, but we also see that the Ai?-distribution 
is rather asymmetric, in contrast to the distribution of 
ln{MAE). The average {\n{M aEq)) decreases as const — 
M/^, with ^ « 29, close to the conductance localization 
length. 

In Fig. H we present our results, i.e. MaE versus M, 
for a repulsive interaction, V = 1.0. The 40 (50,60) sites 
systems have been calculated using m = 190 (375,375) 
states and performing three finite lattice sweeps. The 
accuracy is better than I0~^ for the groundstate energy 
of the 60-sites systems. The phase sensitivity, aE, is 
positive for all samples. The dashed line, obtained by 
fitting (ln(Af Ai?)) as described above, represents —M/^ 
with ^ « 14, about half of the free-fermion value. The 
error in ^ may be about 20%. (We feel that for M = 60, 
the number of samples used, ~ 50, may not be sufficient.) 
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The short-dashed line shows, for comparison, the decay 
of {\n{M aE)) for the non-interacting case. 

In contrast, the phase sensitivity is strongly enhanced 
for an attractive interaction, as can be seen from Fig. ^ 
{V = —1.0). The dashed line, again, is obtained by fit- 
ting (ln(A/A£')) with -M/C, leading to ^ « 100. As an 
example, we have studied a 40-sites system in more de- 
tail, for a repulsive interaction {V = 1.2). Using « 350 
samples, we find a reasonably smooth distribution of 
ln(MAi?)-values, as shown in Fig. ^. These data can 
be fitted well with the corresponding Gaul3 curve, com- 
puted from the average and the variance as given in the 
figure caption; these values were computed from their 
definition. 



CONCLUSIONS 

In conclusion, we have demonstrated that using the 
DMRG-algorithm, it is possible to obtain very accurate 
results for the ground-state properties of interacting one- 
dimensional systems with defects, as e.g. is apparent from 
the comparison of numerical and analytical results for 
the decay of the Friedel oscillations. In particular, sys- 
tem sizes of a few hundred sites are sufficient to obtain 
asymptotic results, except when the decay is very slow, 
i.e. ~ a;~^/^ or slower. Our results for the phase sensi- 
tivity are consistent with what is expected for spinless 
fermions (see ||l2| and reference therein). The attrac- 
tive and the repulsive ground states both contain density 
fiuctuations, but the attractive ground state also contains 
superconducting fluctuations. The latter screen the dis- 
order, leading to an increase of the phase senstivity, i.e. 
a localization length larger than in the non-interacting 
case, for an attractive interaction. On the other hand, 
the phase senstivity is reduced for a repulsive interac- 
tion. As pointed out in |l2), however, these trends are 
just opposite to what should be expected for "realistic" 
models, i.e. models which tend to homogenize the den- 
sity for a repulsive interaction (as e.g. is the case for the 
Hubbard model). We plan to study this question fur- 
ther. We are not aware, however, of any evidence that 
the phase sensitivity can become larger than in the clean, 
non-interacting case. 
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FIG. 1. Decay of the Friedel oscillations induced by weak 
impurities, located summetrically at the ends of the chain 
(ei = — £M = 0.01). The increase for large j arises due to the 
finiteness of the chain. The calculations are performed at half 
filling, N = M/2, keeping m = 120 (200) states per block for 
M — 200 (500). The amplitude of the oscillation vanishes for 



FIG. 2. The exponent 5 vs. interaction, for the same impu- 
rity as in Fig. 11]. The continous line is the asymptotic result 



FIG. 3. Phase sensitivity of the ground state energy vs. 
system size, for a repulsive interaction {V — 1.0, W = 2.0, 
N = M/2). This value of W corresponds to ~ 26 for the 
non-interacting case. The decay length is ^ ~ 14. 

FIG. 4. Phase sensitivity versus system size for an attrac- 
tive interaction {V = -1.0, W = 2.0, = M/2). The decay 
length is 5 100. 

FIG. 5. Distribution of In(MAB) for V = 1.2, M = 40, 
on the basis of 353 samples. We find (In MaE) = —2.54 and 
(7 = 1.5 from our data. The continous line is the correspond- 
ing Gaufi curve. 
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